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Although today there are many observational methods, Type la supernovae (SNIa) is still one of 
the most powerful tools to probe t he mysterious dark en ergy (DE). The most recent SNIa datasets 
are the 307 SNIa " Union" dataset iKowalski et al. 200SH and the 397 SNIa "Constitution" dataset 
[Hicken et al. 200ff ] . In a recent work [Wei 2010( 1 , Wei pointed out that both Union and Constitution 
datasets are in tension with the observations of cosmic microwave background (CMB) and baryon 
acoustic oscillation (BAO), and suggested that two truncated versions of Union and Constitution 
dataset s, namely "UnionT" and "ConstitutionT" , should be used to constrain various DE models. 
But in [Wei 20fo| . only the ACDM model is used to select the outliers from the Union and the 
Constitution dataset. In principle, since different DE models may select different outliers, the 
truncation procedure should be perfor med for ea ch different DE model. In the present work, by 
performing the truncation procedure of [Wei 20ibT] for 10 different mo dels, we d emonstrate that the 
impact of different models is negligible, and the approach adopted in [Wei 2010(| is valid. Moreover, 
by using the 4 SNIa datasets mentioned above, as well as the observations of CMB and BAO, we 
perform best-fit analysis on the 10 models. It is found that: (1) For each DE model, the truncated 
SNIa datasets not only greatly reduce Xmin an d Xmin/dof , but also remove the tension between 
SNIa data and other cosmological observations. (2) The CMB data is very helpful to break the 
degeneracy among different parameters, and plays a very important role in distinguishing different 
DE models. (3) The current observational data are still too limited to distinguish all DE models. 



I. INTRODUCTION 



Observations of Type la supe rnovae (SNIa) [Riess et al. 19981 iPerlmutter et al. 1999L Tonrv et al. 2003 , 
Knop et al. 20031 iRiess et al. 2004al. cosmic microwave backg ro und (CMB) iBennet et al. 20031 ISpergel et al. 20031 



Spergel et al. 20071 iPage et al. 20071 iHinshaw et al. 20071 iKomatsu et al. 20091 and large scale struc- 
ture (LSS) [Tegmark et al. 2004al iTegmark et al. 200411 iTegmark et al. 2006( 1 al1 indicate the existence 
of dark energy (DE) driving the current accelerating expansion of the universe. The most obvi- 
ous theoretical candidate of DE is the cosm ological consta n t A, but it i s plagued with the fine- 
tuning prob le m and the coincidence problem IWeinberg 19891 Weinberg 2000[ ISahni and Starobinskv 2000l 
ICarroll 200il iPeebles and Ratra 20031 iPadmanabhan 20031 ICopeland et al. 20061 There are also 

many dynamical DE models, such a s quintessence IPeebles and Ratra 19881 Ratra and Peebles 19881 



Ratra and Peebles 19881 IZlatev et al. 1999). phantom [CaTdwell 2002[ jcarroll et al. 20031 . fc-essence 



Armendariz- Picon et al. 19991 IChiba et al. 2000l lArmendariz-Picon et al. 200ll. CPL IChevallier and Pola rski 2001, 
iLinder 2003j |, tac hyon IPadmanabhan 20021 iBagla et al. 2003) . hessence IWei et al. 20051 IWei and Cai 2005) . 



Chaplyg in gas |Kamenshchik et al. 200 ill , generali ze d Chaplygin gas Bento et al. 2002 1 , hol ographic 
[Li 20041 iHuang and Li 2005al iHuang and Li 2005bl iLi et al. 2001 \U et al 2009al IZhang 2010) . age- 



graphic IWei and Cai 2008ai IWei and Cai 2008bl . hologra p hic Ricci Gao et al. 20091 . Yang- Mills condensate 
(Zhang et al. 20071 IXia and Zhang 20071 IWang et al. 20081 IWang and Zhang 2008) .etc. Although numerous 
theoretical models have been proposed in the past decade, the nature of DE still remains a mystery. 

In recent years, the numerical study of DE, i.e. utilizing cosmological ob servations to constra i n DE mod- 
els, has become one of t he most active fields in t h e modern cosmology lAlbrecht et al. 2006). Wang 2 000. 
Wang and Mukheriee 2004 I Wang and Mukheriee 20071 IHuterer and Starkman 20031 iHuterer and Coorav 2005 , 



Huang and Gong 2 004. Z hang and Wu 2004 IZhang 20051 IChang et al. 20061 IZhang 20071 IZhang and Wu 2007 
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iGong and Li 2010ai IGong et al. 2010bl ISanchez et al. 20091 iMortonson et al. 20091 IWang et al. 2010} . However, 
these SNIa datasets are not always consist ent with other types of cosmological observations, and are even in tension 
with other SNIa samples. For examples, in [Jassal et al. 20051 1 Jassal et al. 20061 iNesseris and Perivolaropoulos~2 005] . 
the Gold04 dataset was shown to be inconsistent with the SNLS datase t: the SNLS dataset favors the AC DM 
model, while the Gold04 dataset favors the dynamical DE model. In [Ne sseris and Perivolaropoulos 20071 ] . by 
comparing the maximum likelihood fits of the CPL parameters (u>o, wi) given by different SNIa samples, Nesseris 
and Perivolaropoulos found that the Gold06 dataset is also in 2a tension with the SNLS dataset. Moreover, they also 
investigated how to remove this tension. The method is simple. First, they fitted the ACDM model to the whole 182 
SNIa in the Gold06 dataset, and obtained the best-fit parameter value of ACDM model (for SNIa data only). Then, 
they calculated the relative deviation to the best-fit ACDM prediction, \p b s — Pacdm\/ '&obs, for all the 182 data 
points. Here fi b s is the observational value of distance modulus, pacdm is the theoretical value of distance modulus 
given by the best-fit ACDM model, and cr b s is the la error of distance modulus. By searching the SNIa samples 
satisfying \p b s ~ Pacd m \ /&obs > 1-8, they isolated six SNIa that are mostly responsible for the tension. Further, 
by using the random truncation method, they demonstrated that these 6 SNIa are systematically different from the 
Gold06 dataset. 

In a recent work, by comparing the maximum likelihood fits of the CPL parameters (wq, u>i), Wei Wei 2010| pointed 
out that both Union and Constitution dataset are also in tension with the observations of CMB and baryon acoustic 
oscill ation (BAO). Moreover, he also in vestigated how to remove these tensions. By using the method of truncation 
of [Nesseris and Perivolaropoulos 20071 ] . Wei found out the main sources that are responsible for the tensions: for the 
Union set, there are 21 SNIa differing from the best-fit ACDM prediction beyond 1.9cr (i.e. \p bs~ HACDM\/<Jobs > 1-9); 
and for the Constitution set, there are 34 SNIa differing from the best-fit ACDM prediction beyond 1.9a. The specific 
limit of truncation (i.e. 1.9er) is chosen based on two considerations: first, the tension between SNIa samples and other 
obser vations ca n be completely removed; second, the number of usable SNIa can be preserved as much as possible 
(see [Wei 2010| for details). By subtracting these outliers from the Union dataset and the Constitution dataset, 
respectively, two new SNIa datasets, "UnionT" and "ConstitutionT" ("T" stands for "truncated"), were obtained. 
Further, Wei argued that the UnionT and the ConstitutionT datasets are f ully consis tent with the other cosmological 
observations, and should be used to constrain various DE models. But in [Wei 2010} . only the ACDM model is used 
to select the outliers from the Union and the Constitutio n dataset. Since different DE models may select different 
outliers, one may doubt whether the approach adopted in [Wei 20foj is valid. In principle, the truncation procedure 
should be performed for each different DE model. Only if the impact of different models is negligible, one can conclude 
that the conclusion of Wei is correct. So in the present work, we shall consider 10 different models. The truncation 
procedure will be performed for all these 10 models, and the corresponding cosmological consequences will be explored. 

This paper is organized as follows: In Section 2, we briefly describe 10 theoretical models considered in this work. 
In Section 3, we present the method of data an alysis, as well as the SNIa datasets we used in this paper. By 
pe rforming th e truncation procedure of [Wei 2010l j for 10 different models, we demonstrate that the approach adopted 



Wei 2010] is valid. In Section 4, we show the data fitting results of 10 models, and present the corresponding 
conclusions. Section 5 is a short summary. In this work, we assume today's scale factor do = 1, so the redshift z 
satisfies z = a -1 — 1; the subscript "0" always indicates the present value of the corresponding quantity, and the unit 
with c = h — 1 is used. 



II. MODELS 



For a spatially flat (the assumption of flatness is motivated by the inflation scenario) Friedmann- Robertson- Walker 
(FRW) universe with matter component p m and DE component pde, the Friedmann equation reads 



Ml 2 Pl H 2 = Pm + p de , 



(1) 



3 



where Mpi = l/vo7rG is the reduced Planck mass, and H = a/a is the Hubble parameter. Using this formula, one 
can easily get 

E(z) ee H(z)/H Q = [n m0 (l + zf + (1 - n mQ )f(z)] 1/2 , (2) 

where Hq is the Hubble constant, f2 m o is the present fractional matter density, and key function f(z) = pde{%) / Pde(O) 
is given by the specific DE model. Equivalently, we have 

where fide = ^f- = zm^h 12 * s ^ ne fractional DE density. Clearly, different DE model will give different E(z). In the 
following, we shall briefly describe the models considered in this work. 

A. Single-parameter models 

(1) The ACDM model: The DE density is always a constant, i.e. the equation of state (EOS) w = —1, so 



e(z) = y/n m o{i + z) 3 + (i-n m0 ). (4) 

(2) The Dvali -Gabadadze-Porrati (DGP) model: For this modified gravity model, the form of E(z) is 
jDvali et al. 2000j ] 



e{ z ) = y/n m0 (i + z) 3 + n rc + v^T, (5) 

where tt rc is a constant, satisfies 

^=(^^) 2 - (6) 

(3) The agegraph ic dark energy (ADE): The DE density is characterized by the conformal age r\ of the universe 
(Wei and Cai 2008bj . 

2 r 2 -2 f dt f da 

p de = 3n M Pl rj , and r] = — = (7) 

where n is a positive constant. The equation of motion for fide is given by |Wei and Cai 2008bl | 

(8) 



— 7— — — ^de (1 — ^de) 

dz 



3(1 + z)- 1 - -v 7 ^ 
n 



As in [Wei and Cai 2008bj . we choose the initial condition, Qde(zini) = n 2 (l + z ini )~ 2 /4:, at z ini — 2000, then Eq. ([8]) 
can be numerically solved. Substituting the results of Eq. ((8]) into Eq. ((3]), the key function E(z) can be obtained. 
Notice that once n is given, ft m o = 1 — £!de(0) can be naturally obtained by solving Eq. ((8]), so the ADE model is a 
single-parameter model. 

B. Two-parameter models 

(4) The XCDM model: DE has a constant EOS w, then 



E(z) = ^Jn m0 (l + zf + (1 - n m0 )(l + z)3U+«0. (9) 
(5) The Chaplygin gas (CG) model: The EOS of DE has the form [Kamenshchik et al. 200f| 



Pde = -—, (10) 

Pde 



where A is a positive constant, and pde is the pressure of DE. From this assumption, one can get 
Kamenshchik et al. 2001 



e(z) = sj n m0 (i + zf + (i - n m0 )^A c + (i - a c )(i TW, 



(ii) 



here A c = A/pd e (0) is also a positive constant. 

(6) The ho lographi c dark energy (HDE) model: The DE density is characterized by the future event horizon L of 
the universe Li 20041 . 



f°° dt' 

Pde = 3c 2 Mp,L~ 2 , and L = a — = a 

Jt a 



da' 



where c is a positive constant. The equation of motion for Q^e is given by [Li 2004 1 

dQ,de x _1 „ ,.. „ s /„ 2 



= -(l + z^n^i-nde) (1 + 



(12) 



(13) 



Solving Eq. (fl3|) numerically and substituting the corresponding results into Eq. (j^J, £'(z) can be obtained. 

(7) The holog raphic Ricci dark energy (RDE) model: The DE density is characterized by the Ricci scalar 1Z 



Gao et al. 200 



Pde 



TZ, and K=-g(h + 2H 2 ^ , 



where a is a positive constant. The form of E(z) in this case is Gao et al. 2009| 

E{z) = 



^(1 + zf + (1 - ^)(1 + z)( 4 -^. 



2- a 



2-a' 



(14) 



(15) 



C. Three- parameter models 

(8) The linear parameterization (LP) model: The EOS of DE is parameterized as 

w — w + WiZ, 

where wq and uii are constants. Then 



E{z) = yt(l + ^) 3 + (1 - fi m0 )(l + z)3(i+-o-- 1 ) exp(3 Wl 2). 



(16) 



(17) 



(9) The C hevallier-Polarski-Linder (CPL) model: The EOS of DE is parameterized as [Chevallier and Polarski 20011 
iLinder 20031 ] 



W = Wq + Wi 



l + z' 



then 



E(z) = yd(l + zf + (1 - fi m0 )(l + z )3(i +u , +«-i) exp ■ 



(18) 



(19) 



(10) The generalized Chaplygin gas (GCG) model: The EOS of DE has the following form [Bento et al. 2002j | 

A 



Pde 



Pde 



(20) 



where a is also a positive constant and a = 1 corresponds to the CG model. One can get [Bento et al. 2002 1 

E(z) = vWo(l + zf + (1 - fi m0 ) (As + (1 - A s )(l + z)3(i+«)) 1/1+a . (21) 



here A s = A/p d + e a {0) 
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III. METHODOLOGY 
A. Data analysis 

In this work we adopt \ 2 statistics. For a physical quantity £ with experimentally measured value £ bs, standard 
deviation er^, and theoretically predicted value £th, the \ 2 value is given by 

2 _ (Ct/l _ Cobs) fr,r,\ 

Xc -5 ■ v zz ) 



The total \ 2 is the sum of all xf s , i- e - 



x 2 



= £4 ( 23 ) 



First, we consider the SNIa data that are given in terms of the distance modulus n b s (zi). The theoretical distance 
modulus is defined as 

Hth(zi) = 51og 10 D L (zi) + fi , (24) 

where fio = 42.38 — 51og 10 h with h the Hubble constant Hq in units of 100 km/s/Mpc, and in a flat universe the 
Hubble- free luminosity distance Dl = H cIl (<1l denotes the physical luminosity distance) is 

W'l'^'l'w (25) 

where 9 denotes the model parameters. The x 2 f° r the SNIa data is 

Xsn(9) = 2^ 2 > V 26 ) 

where ix bs{zi) and <7j are the observed value and the corresponding lcr error of distance modulus for each supernova, 
respectively. Paramete r is a nuisance parameter but it is independent of the data and the dataset. Following 



Perivolaropoulos 2005j , the minimization with respect to fj,o can be made trivially by expanding the \ 2 °f Eq. (|26|) 
with respect to \i§ as 



where 



X 2 SN (0) = A(9) - 2 l x Q B(0) + fiC, (27) 

A ^ _ [v-obsjzj) - nth{zi; Mo = 0,9)] 2 ^ 

B ( ) = ~ /itfe ( z " Mo I (29) 

C = Y,^2- (30) 
Evidently, Eq. (f26)) has a minimum for ^io = B/C at 

^(0)=A(0)-^. (31) 

Since x'sNmin = Xsn min> ms tead of minimizing Xsn we w ^ minimize Xsn which is independent of the nuisance 
parameter (Iq. 
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Next, we consider co nstraints from the CMB and the LSS observat ions. For the CMB data, we use the CMB shift 
parameter R, given by Bond et al. 19971 IWang and Mukheriee 2006j 



R = n 



1/2 
raO 



dz' 



(32) 



1091.3 iKomatsu et al. 20101 



The m easured value of R has been updated 
2010} . For the LSS data, we use the BAO 



where the redshift of recombination z rec 

to be R obs = 1.725±0.018 from the WMAP7 observations [Komatsu et al. 
distance meas urements obtained at z = 0.2 and z = 0.35 from the joint analysis of the 2dFGRS a nd SPSS data 
(Percival 2009j . The BAO distance ratio D v (z = 0.35) /D v (z = 0.20) = 1.736 ± 0.065 was shown in jPercival 2009| 
to be a relatively model independent quantity. Here Dy(z) is defined as 



D v (z B Ao) = 



zbao 
H(zbao) 



dz 
Wz) 



1/3 



The total % 2 is given by 



~2 2 2 

XSN + XCMB + XBAO > 



where Xsn i s given by Eq. pip , and the latter two terms are defined as 

(R- 1.725) 2 



XCMB 



0.018 2 



and 



2 _ (L>y(0.35)/L>y(0.20) - 1.736) 2 
Xbao — 



0.065" 



(33) 



(34) 



(35) 



(36) 



The model parameters yielding a minimal x 2 is favored by the observations. 

To compare different models, a statistical variable must be chosen. The Xmin 1S the simplest one, but it has difficulty 
to compare different models with different number of parameters. In this work, we will use Xmin/dof as a model 
selection criterion, where dof is the degree of freedom defined as 



dof = N- k, 



(37) 



here N is the number of data, and k is the number of free parameters. This model selection criterion has been widely 
used in the literature. 



B. SNIa Datasets 



Although today there are many observational methods, SNIa is still one of the most powerful tools to probe the 
mysterious DE. In the past decade, many SNIa datasets have been released, w hile the number and quality of SNIa 
have continually increased. In 2008, the Union dataset [Kowalski et al. 2008| was released. It includes the large 
samples of SNIa from the HST, SNLS and ESSENCE, and contai ns 307 samples: 250 high redshift SNIa (z > 0.2) 
and 57 low redshift SNIa (z < 0.2). In 2009, the Constitution set jHicken et al. 2009j was released. It adds 90 low 
redshift samples. Th ese two SN Ia datasets have been widely used in the literature. 

In a recent work Wei 2O10j . by comparing the maximum likelihood fits of the CPL parameters (w , Wi), Wei 



pointed out that both Union and Constitution datasets are in tension not only with the observations of CMB and 
BAO, but also with the other SNIa datasets. Moreover, he also investigated how to remove these tensions. The 
method of Wei is as follows. First, he fitted the ACDM model to all the data points in the specific SNIa dataset, 
and obtained the best-fit parameter value of ACDM model (for SNIa data only). Then, he calculated the relative 
deviation to the best-fit ACDM prediction, \fJ, bs ~ ^ACDm\/ '&obsi for all the SNIa data points. By searching the SNIa 
samples satisfying \p, bs — ^kcdm\/ &obs > 1-9, Wei get the main sources that are responsible for the tensions. For the 
Union set, there are 21 SNIa differing from the best-fit ACDM prediction beyond 1.9er, which is called "UnionOut" 
subset. 

• UnionOut subset (21 SNIa): 

1992bs, 1995ac, 1999bm, 1997o, 2001hu, 1998ba, 04Pat, 05Red, 2002hr, 03D4au, 
04D3cp, 03Dlfc, 03D4dy, 03Dlco, bOlO, d033, g050, g055, k430, ml38, m226 
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TABLE I: The "UnionOut" samples that satisfying \fi bs — ^modei\/ &obs > 1.9, for each model. Here "+" denotes adding a data 
point, while "-" denotes reducing a data point. 



ivioQei 


Numbers of "UnionOut" samples for each model 


Differences with the ACDM UnionOut subset 


ACDM 


21 




None 




DGP 


21 




None 




A DP 

JrVUAlj 


Z 1 




None 




XCDM 


22 


+1999gd 


+f076 


-05Red 


CG 


22 


+1999gd 


+f076 


-05Red 


HDE 


22 


+ 1999gd 


+f076 


-05Red 


RDE 


22 


+1999gd 


+f076 


-05Red 


LP 


19 


+1999gd +f076 -1999bm -05Red -03Dlco -k430 


CPL 


20 


+1999gd +f076 - 


1999bm 


-05Red -03Dlco 


GCG 


22 


+1999gd 


+f076 


-05Red 



For the Constitution set, There are 34 SNIa differing from the best-fit ACDM prediction beyond 1.9a, which is called 
"ConstitutionOut" subset. 

• ConstitutionOut subset (34 SNIa): 

1992bs, 1992bp, 1995ac, 1999bm, 1996t, 1997o, 1995aq, 2001hu, 1998ba, 04Pat, 05Red, 
2002hr, 03D4au, 04D3gt, 04D3cp, 03D4at, 03Dlfc, 04D3co, 03D4dy, 04D3oe, 04Dlak, 
03Dlco, bOlO, d033, f076, g050, k430, ml38, m226, snOlcp, sn02hd, sn03ic, sn07ca, sn07R 

It should be mentioned that 1.9ct is chosen based on two considerations: first, the tension between SNIa samples 
and other observations can be completely removed; second, the number of usable SNIa can be preserved as much as 
possible. After taking into account these two factors, Wei found that 1.9cr is most appropriate for the truncation 
procedure. By subtracting these two subsets from the Union dataset and the Constitution dataset, respectively, two 
new SNIa samples, "UnionT" and "ConstitutionT" , were obtained. It is clear that UnionT has 286 SNI a samples , 
and ConstitutionT has 363 SNIa samples. Analyzing these two truncated datasets with CPL model, Wei Wei 2O10j 
argued t hat they ar e fully consistent with the other cosmological observations. 

But in [Wei 20 10]. only the ACDM model is used to select the outliers from the Union and the Constitution dataset. 
Since different DE models may select different o utliers (i.e . different "UnionOut" and "ConstitutionOut" Samples), 
one may doubt whether the approach adopted in Wei 2010l | is valid. In principle, the truncation procedure should be 



performed for each different DE model. Only if the impact of different models is negligible, one can conclude that the 
conclusion of Wei is correct. So in this work, we perform the truncation procedure for 10 different DE models. As in 
Wei 2010j . we choose 1.9tr as the selection criterion. The results are shown in table Hand table HH From these two 



tables, it is seen that the difference among the outlier samples given b y various m odels are very small, and the impact 
of different models is negligible. Therefore, the approach adopted in [Wei 2010j is valid. Since the difference among 
the outlier samples given by various mode ls are very small, for simplicity, in this work we just adopt the UnionT and 
the ConstitutionT SNIa samples given by Wei 2010] . To explore the corresponding cosmological consequences, in the 



following we shall study 10 models by using 4 SNIa datasets: Union, UnionT, Constitution and ConstitutionT. 



IV. RESULTS AND CONCLUSIONS 



We will present the results of data fitting in this section. Using the Union, the UnionT, the Constitution and the 
ConstitutionT dataset, respectively, we list the xLin an d the Xmin/dof for those 10 models, in table IIII1 table IIV1 
table[V] and table I VII Here "SNIa" means only SNIa data is used in the analysis, "SNIa+BAO" means both SNIa data 
and BAO data are used, "SNIa+CMB" means both SNIa data and CMN data are used, and "SNIa+BAO+CMB" 
means all these three types of observational data are taken into account. The main conclusions are summarized as 
follows. 

(1) For each DE model, the truncated SNIa datasets not only greatly reduce Xmin an d Xmin/d°fi but also remove 
the tension between SNIa data and other cosmological observations. 
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TABLE II: The "ConstitutionOut" samples that satisfying \fi bs — jJ-modei\/oobs > 1.9, for each model. Here "+" denotes adding 
a data point, while "-" denotes reducing a data point. 



ivioaei 


Numbers of 


"ConstitutionOut" samples for each model Difference with the ACDM ConstitutionOut subset 


ACDM 




34 






None 


DGP 




34 






None 






33 






-sn07R 


XCDM 




34 






None 


CG 




34 






+1997aj -1992bp 


HDE 




34 






None 


RDE 




34 






None 


LP 




32 




+1997aj 


-1992bp -1999bm -03D4at 


CPL 




33 




+ 1997aj +sn02bf -1992bp -1999bm -sn07R 


GCG 




34 






+1997aj -1992bp 


TABLE III: The xlir, 


and the Xmin/dof ( 


in the Parentheses) for the 10 models, where the Union dataset is used. 




Model 


SNIa 


SNIa+BAO 


SNIa+CMB 


SNIa+BAO+CMB 




ACDM 


311.936 (1.019) 


313.205 (1.020) 


312.424 (1.018) 


313.594 (1.018) 




DGP 


313.026 (1.023) 


314.319 (1.024) 


339.280 (1.105) 


341.584 (1.109) 




ADE 


313.536 (1.025) 


314.838 (1.026) 


327.216 (1.066) 


329.252 (1.069) 




XCDM 


310.682 (1.019) 


311.966 (1.020) 


312.224 (1.020) 


313.456 (1.021) 




CG 


310.434 (1.018) 


311.628 (1.018) 


311.921 (1.019) 


313.193 (1.020) 




HDE 


310.827 (1.019) 


312.149 (1.020) 


311.218 (1.017) 


312.481 (1.018) 




RDE 


310.682 (1.019) 


311.966 (1.020) 


312.323 (1.021) 


313.801 (1.022) 




LP 


309.984 (1.020) 


310.744 (1.019) 


310.691 (1.019) 


311.978 (1.020) 




CPL 


310.091 (1.020) 


310.896 (1.019) 


310.906 (1.019) 


312.258 (1.020) 




GCG 


310.405 (1.021) 


311.401 (1.021) 


311.925 (1.023) 


313.325 (1.024) 


TABLE IV: The X L„ 


and the Xmin/dof (in the Parentheses) for the 10 models, where the UnionT dataset is used. 




Model 


SNIa 


SNIa+BAO 


SNIa+CMB 


SNIa+BAO+CMB 




ACDM 


204.568 (0.718) 


205.945 (0.720) 


205.575 (0.719) 


206.794 (0.721) 




DGP 


205.234 (0.720) 


206.634 (0.722) 


226.987 (0.794) 


229.388 (0.799) 




ADE 


205.560 (0.721) 


206.963 (0.724) 


216.201 (0.756) 


218.305 (0.761) 




XCDM 


204.058 (0.719) 


205.449 (0.721) 


204.855 (0.719) 


206.207 (0.721) 




CG 


204.050 (0.718) 


205.365 (0.721) 


204.552 (0.718) 


205.930 (0.720) 




HDE 


204.070 (0.719) 


205.496 (0.721) 


204.228 (0.717) 


205.614 (0.719) 




RDE 


204.058 (0.719) 


205.449 (0.721) 


205.875 (0.722) 


207.498 (0.726) 




LP 


204.055 (0.721) 


205.347 (0.723) 


204.429 (0.720) 


205.739 (0.722) 




CPL 


204.057 (0.721) 


205.391 (0.723) 


204.063 (0.719) 


205.517 (0.721) 




GCG 


204.053 (0.721) 


205.366 (0.723) 


204.545 (0.720) 


206.090 (0.723) 
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TABLE V: The Xmin an d the Xmin/dof (in the Parentheses) for the 10 models, where the Constitution dataset is used. 



Model 


SNIa 


SNIa+BAO 


SNIa+CMB 


SNIa+BAO+CMB 


ACDM 


4o0.oU4 


(I. Kb) 


466.902 (1.176) 


466.316 (1.175) 


467.525 (1.175) 


DGP 


400.1zz 


( 1 1 r 7 r 7\ 
(1.1 / I ) 


467.433 (1.177) 


498.264 (1.255) 


500.368 (1.257) 


ADE 


A(\P> 97^ 


( 1 1 77~\ 
(1.1/ / ) 


467.580 (1.178) 


483.675 (1.218) 


485.585 (1.220) 


XCDM 


465.602 


(1.179) 


466.901 (1.179) 


465.657 (1.176) 


466.947 (1.176) 


CG 


465.293 


(1.178) 


466.564 (1.178) 


465.591 (1.176) 


466.891 (1.176) 


HDE 


465.769 


(1.179) 


467.080 (1.179) 


466.030 (1.177) 


467.384 (1.177) 


RDE 


465.602 


(1.179) 


466.901 (1.179) 


472.841 (1.194) 


474.466 (1.195) 


LP 


461.071 


(1.170) 


461.570 (1.169) 


465.610 (1.179) 


466.910 (1.179) 


CPL 


461.526 


(1.171) 


462.112 (1.170) 


465.636 (1.179) 


466.902 (1.179) 


GCG 


465.080 


(1.180) 


466.360 (1.181) 


465.920 (1.180) 


466.895 (1.179) 


The Xmin 


and the Xmin/ do f (in 


the Parentheses) for the 10 models, where the ConstitutionT dat< 


Model 


SNIa 


SNIa+BAO 


SNIa+CMB 


SNIa+BAO+CMB 


ACDM 


269.081 


(0.743) 


270.610 (0.745) 


271.423 (0.748) 


272.764 (0.749) 


DGP 


269.443 


(0.744) 


270.979 (0.747) 


292.067 (0.805) 


294.360 (0.809) 


ADE 


269.613 


(0.745) 


271.139 (0.747) 


280.090 (0.772) 


282.147 (0.775) 


XCDM 


269.066 


(0.745) 


270.599 (0.748) 


269.245 (0.744) 


270.754 (0.746) 


CG 


268.992 


(0.745) 


270.502 (0.747) 


269.066 (0.743) 


270.594 (0.745) 


HDE 


269.113 


(0.745) 


270.664 (0.748) 


269.130 (0.743) 


270.698 (0.746) 


RDE 


269.066 


(0.745) 


270.599 (0.748) 


273.305 (0.755) 


275.180 (0.758) 


LP 


268.812 


(0.747) 


270.021 (0.748) 


269.068 (0.745) 


270.603 (0.748) 


CPL 


268.900 


(0.747) 


270.147 (0.748) 


269.122 (0.745) 


270.676 (0.748) 


GCG 


268.980 


(0.747) 


270.473 (0.749) 


269.069 (0.745) 


270.599 (0.748) 



Since it is too prolix to describe the results of all 10 models, in the following we will just discuss the HDE 
model as an example. From table IIIII and table IIV1 we can see that the combined Union+BAO+CMB data gives 
a Xmin = 312.481 and a x min /dof = 1-018, while the combined UnionT+BAO+CMB data gives a x min = 205.614 
and a Xmin/dof — 0.719. This means that by dropping the 21 outliers, The Xmin an( l l ne Xmin/ do f °f l ne HDE 
model can reduce 106.867 and 0.299, respectively. From table [V] and table IVI1 similar results are obtained. The 
combined Constitution+BAO+CMB data gives a Xmin — 467.384 and a Xmin/dof = 1-177, while the combined 
ConstitutionT+BAO+CMB data gives a Xmin = 270.698 and a Xmin/ do f — 0.746. This means that by dropping 
the 34 outliers, The x m in an( l t ne Xmin/dof of the HDE model can reduce 196.686 and 0.431, respectively. It 
should be pointed out that for all these 10 models, the UnionT and the ConstitutionT dataset can greatly reduce the 
corresponding Xmin an d Xmin/dof- Moreover, the decreasing margin of X m in an( l Xmin/dof for these 10 DE models 
are almost same. 



To further verify this conclusion, we also adopt the method of random truncation used 



Nesseris and Perivolaropoulos 20071 I Wei 201(ij . The method is very simple. First, we random select the out- 
liers from the full Union dataset and the full Constitution dataset, respectively. Notice that the numbers of data 
points in these random "UnionOut" subsets are same as that in the original UnionOut samples, while the numbers of 
da ta points in these random "ConstitutionOut" su bsets are same as that in the original ConstitutionOut samples. As 
in Nesseris and Perivolaropoulos 20071 IWei 20l6| , 500 random UnionOut subsets and 500 random ConstitutionOut 



subsets are selected. By subtracting these outliers from the Union dataset and the Constitution dataset, respectively, 
500 random UnionT subsets and 500 random ConstitutionT subsets are obtained. Then, we perform best-fit analysis 
on the HDE model by using these random truncated SNIa subsets. Using the SNIa data only, we get the Mean 
of Xmin l° r 500 random UnionOut subsets and 500 random ConstitutionOut subsets, respectively. Notice that the 
full Union set gives a Xmin — 310.827 and the full Constitution set gives a Xmin = 465.769, the differences between 
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TABLE VII: A comparison for the original truncated SNIa subset and the random truncated SNIa subsets. Notice that the full 
Union set gives a Xmin = 310.827 and the full Constitution set gives a Xmin = 465.769. Here AN is the number of data points 
in outliers, x 2 i s the Mean of Xmin f° r those random truncated SNIa subsets, and Ax 2 is the differences between the Xmin °f 
the full SNIa dataset and the Xmin (° r the Mean) of the random truncated SNIa subsets. 



The subset of outliers 


AN 


2 — 2 

Xmin Or X 


Ax 2 


A X 2 /AN 


The original UnionOut subset 


21 


204.070 


106.757 


5.084 


500 random UnionOut subsets 


21 (Mean) 


289.617 (Mean) 


21.210 (Mean) 


1.010 (Mean) 


The original ConstitutionOut subset 


34 


269.113 


199.656 


5.784 


500 random ConstitutionOut subsets 


34 (Mean) 


425.763 (Mean) 


40.006 (Mean) 


1.177 (Mean) 



the Xmin °f the full SNIa dataset and the Mean of Xmin °f these random truncated SNIa subsets are also obtained. 
Next, we make a comparison for the original truncated SNIa subset and the random truncated SNIa subsets. The 
results are shown in table I VIII It is found that the Xmin °f the HDE model can reduce 5 by dropping 1 data point in 
the original truncated SNIa subset, but can only reduce 1 by dropping 1 data point in the random truncated SNIa 
subset. Therefore, the original UnionOut and ConstitutionOut subsets are systematically different from the full 
Union and Constitution datasets, and the original truncated supernova datasets provide a significantly better model 
fittin g than the full SNIa datasets. 

In [Wei 20l"o| , by performing the best- fit analysis on the CPL model, Wei argued that the UnionT and the Con- 
stitutionT datasets are fully consistent with the other cosmological observations. To verify this conclusion, more 
theoretical models should be taken into account. Here we study the HDE model as an example. In Fig|TJ we plot 
the la and the 2a confidence level (CL) contours for the HDE model, where the Union and the UnionT datasets are 
used, respectively. From this figure, we find that the best-fit point for the Union data is outside the 2a confidence 
region given by the combined Union+BAO+CMB data, while the best-fit point for the UnionT data is inside the 
la confidence region given by the combined UnionT+BAO+CMB data. This means that the UnionT dataset is 
very useful to remove the tension between SNIa data and other cosmological observations. In addition, we also use 
the Constitution and the ConstitutionT datasets to plot the CL contours for the HDE model in FigO It is found 
that the best-fit point for the Constitution data is outside the 2a confidence region given by the combined Consti- 
tution+BAO+CMB data, while the best-fit point for the ConstitutionT data is very close to the best-fit point for 
the combined ConstitutionT +BAO+CMB data. This means that the ConstitutionT dataset is also very helpful to 
remove the tension. Therefore, we conclude that the truncated SNIa datasets can remove the tension between SNIa 
data and other cosmological observations. 

(2) The CMB data is very helpful to break the degeneracy among different parameters, and plays a very important 
role in distinguishing different DE models. 

As an example, we will compare the HDE model with the RDE model. By using the Constitution SNIa data alone, 
the CMB data alone, and the combined Constitution+CMB+BAO data, respectively, we plot the la and the 2a CL 
contours for the HDE and the RDE model in FigJ3] From this figure, we find that the shapes of CL contours given by 
the CMB data alone are quite different from that given by the SNIa data alone, and the CMB data is very helpful to 
break the degeneracy among different parameters. As shown in the left panel, the la CL contour of the HDE model 
given by the CMB data intersects to the la CL contour of the HDE model given by the SNIa data; while in the 
right panel, the la CL contour of the RDE model given by the CMB data does not intersect to the la CL contour of 
the RDE model given by the SNIa data. This fact explains why the HDE model performs much better in fitting the 
combined Constitution+BAO+CMB data than the RDE model. Besides, we also check the effect of the BAO data, 
and find that the la confidence region of the HDE model given by the BAO data alone can completely cover that 
given by the SNIa data. So the constraint given by the BAO data is quite weaker than that given by the SNIa and 
the BAO data. This conclusion can be further verified by using table |V] As seen in table |V] without adopting the 
CMB data, it is very difficult to distinguish the HDE Model from the RDE model. After adding the CMB data, it 
is seen that the x mm of the HDE Model given by the SNIa+CMB data is 6.811 smaller than the x min of the RDE 
Model given by the SNIa+CMB data, while the x mm of the HDE Model given by the combined SNIa+CMB+BAO 
data is 7.082 smaller than the x min of the RDE Model given by the combined SNIa+CMB+BAO data. 

In addition, by using the ConstitutionT SNIa data alone, the CMB data alone, and the combined Constitu- 
tionT+CMB+BAO data, respectively, we also plot the la and the 2a CL contours for the HDE and the RDE model 
in FigJU Again, we see that the la CL contour of the HDE model given by the CMB data intersects to the la 
CL contour of the HDE model given by the SNIa data, while the la CL contour of the RDE model given by the 
CMB data does not intersect to the la CL contour of the RDE model given by the SNIa data. This fact explains 
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FIG. 1: The la and the 2a CL contours for the HDE model. The left panel is plotted by using the Union dataset, while 
the right panel is plotted by using the UnionT dataset. For both these two panels, the blue dashed lines correspond to 
the constraints given by the SNIa data only, and the red solid lines correspond to the constraints given by the combined 
SNIa+BAO+CMB data. Moreover, we also plot the best-fit point for the SNIa data (blue point) and the best-fit point for the 
combined SNIa+BAO+CMB data (red star). It should be pointed out that, in the left panel, the best-fit point for the SNIa 
data is outside the 2a confidence region given by the combined SNIa+BAO+CMB data, while in the right panel, the best-fit 
point for the SNIa data is inside the la confidence region given by the combined SNIa+BAO+CMB data. This means that 
the UnionT dataset is very useful to remove the tension between SNIa data and other cosmological observations. 




FIG. 2: The same as in Fig[TJ except for the cases of the Constitution and the ConstitutionT datasets. It should be men- 
tioned that, in the left panel, the best-fit point for the SNIa only is outside the 2a confidence region given by the combined 
SNIa+BAO+CMB data, while in the right panel, the best-fit point for the SNIa only is very close to the best-fit point for the 
combined SNIa+BAO+CMB data. This means that the ConstitutionT dataset is very helpful to remove the tension between 
SNIa data and other cosmological observations. 



why the HDE model performs much better in fitting the combined ConstitutionT+BAO+CMB data than the RDE 
model. Therefore, the CMB data is very helpful to break the degeneracy among different parameters, and plays a 
very important role in distinguishing different DE models. 

(3) The current observational data are still too limited to distinguish all DE models. 

First, we discuss three kinds of two-parameter models: the XCDM model, the CG model, and the HDE model. 
As seen in table El the Xmin °f these 3 models given by the combined Constitution+BAO+CMB data are 466.947, 
466.891, and 467.384, respectively. That is to say, after taking into account the CMB data, the differences of the Xmin 
of these 3 models are still smaller than 1. Notice that this result also holds true in table UTTl table ITVl and table fVTl 
Therefore, one cannot judge which DE model is better. 

Then, we discuss three kinds of three-parameter models: the LP model, the CPL model, and the GCG model. 
As seen in table [V] the Xmin °f these 3 models given by the combined Constitution+BAO+CMB data are 466.910, 
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FIG. 3: The la and the 2a CL contours for the HDE and the RDE model. The Constitution dataset is used to plot this 
figure. The left panel is plotted by using the HDE model, while the right panel is plotted by using the RDE model. For 
both these two panels, The blue dashed lines correspond to the constraints given by the SNIa data only, the black dotted 
lines correspond to the constraints given by the CMB data only, the red solid lines correspond to the constraints given by the 
combined SNIa+BAO+CMB data, and the red stars denote the best-fit point for the combined SNIa+BAO+CMB data. Since 
the shapes of CL contours given by the CMB data alone are quite different from that given by the SNIa data alone, the CMB 
data is very helpful to break the degeneracy among different parameters. As shown in the left panel, the la CL contour of 
the HDE model given by the CMB data intersects to the la CL contour of the HDE model given by the SNIa data; while in 
the right panel, the la CL contour of the RDE model given by the CMB data does not intersect to the la CL contour of the 
RDE model given by the SNIa data. This fact explains why the HDE model performs much better in fitting the combined 
Constitution+BAO+CMB data than the RDE model. Besides, we also check the effect of the BAO data, and find that the 
la confidence region of the HDE model given by the BAO data alone can completely cover that given by the SNIa data. For 
simplicity, the CL contours given by the BAO data are not plotted in this figure. 
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FIG. 4: The same as in FigO except for the case of the ConstitutionT dataset. Notice that the shapes of CL contours given by 
the CMB data alone are quite different from that given by the SNIa data alone. As shown in the left panel, the la CL contour 
of the HDE model given by the CMB data intersects to the la CL contour of the HDE model given by the SNIa data; while 
in the right panel, the la CL contour of the RDE model given by the CMB data does not intersect to the la CL contour of 
the RDE model given by the SNIa data. This fact explains why the HDE model performs much better in fitting the combined 
ConstitutionT+BAO+CMB data than the RDE model. 



466.902, 466.895, respectively, and the differences of the Xmin °f these 3 models are even smaller than 0.1. Since 
this result also holds true in table IIII1 table IIV1 and table I VII it is also very difficult to distinguish these 3 models. 
Therefore, to distinguish DE models better, more high-quality observational data are needed. 
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V. SUMMARY 



In this work, by perfor ming the truncation procedure of [Wei 2010} for 10 different models, we demonstrate that 
the approach adopted in [Wei 20l61 ] is valid. Moreover, by using the 4 SNIa datasets mentioned above, as well as the 
observations of CMB and BAO, we perform best-fit analysis on the 10 DE models. It is found that: (1) For each 
DE model, the truncated SNIa datasets not only greatly reduce Xmin an d Xmin/^ /' but also remove the tension 
between SNIa data and other cosmological observations. (2) The CMB data is very helpful to break the degeneracy 
among different parameters, and plays a very important role in distinguishing different DE models. (3) The current 
observational data are still too limited to distinguish all DE models. These results provide a further support for Wei's 
work, and indicate that the existence of biasing systematic errors in SNIa data should be taken into account seriously. 
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